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ABSTRACT 

A survey of solar wind ion velocity distributions and derived parameters 
(temperature, ion differential speed, heat flux, adiabatic invariants) is pre- 
sented with emphasis on the heliocentric distance range between 0.3 and 1 AU 
traversed by the Helios solar probe. We discuss the radial evolution of non-ther- 
mal features (e.g., T x > T „ , Av^ = V^), which are observed to be most pro- 
nounced at perihelion?WithiR the - ^framework of quasilinear plasma theory, wave- 
particle interactions that may shape the ion distributions are considered. Some 
results of a self-consistent model calculation are presented accounting for ion 
acceleration and heating by resonant momentum and energy exchange with ion-cy- 
clotron and magnetosonic waves propagating away from the sun along the inter- 
planetary magnetic field. Another tentative explanation for the occurrence of 
large perpendicular proton temperatures is offered in terms of heating by Landau 
damping of lower hybrid waves. 


1. Introduction 

The purpose of this paper is to review some experimental and theoretical 
work that has recently been done on the problem of the radial evolution of solar 
wind ion velocity distributions. We shall concentrate on Helios observations which 
complement previous in situ measurements at 1 AU and at larger heliocentric dis- 
tances with very detailed ion measurements in the inner heliosphere between 0.3 
and 1 AU. For this radial distance range average gradients of fluid parameters 
are now available that put further contraints on theoretical modeling of the solar 
wind expansion. Our results may also help to better understand the evolution of 
ion internal energy and non-thermal characteristics of the distributions. To date, 
the most detailed reviews on this subject have been published by Feldman et al. 
(1974) on observations of interpenetrating solar wind streams and non-thermal 
distributions, by Hollweg (1974) on waves and instabilities, and by Feldman et al. 
(1979a, b) , and Schwartz (1980) on kinetic processes and microinstabilities, re- 
spectively, that affect and shape ion distribution functions. Associated plasma 
waves have been reviewed by Gurnett (1981) and Scarf et al. (1981). This brief 
paper does not aim at a complete treatment of all relevant issues but is biased 
towards the author’s preferences. We apologize at the outset for inevitable omis- 
sions of references. An exhaustive discussion of all topics is beyond the scope 
of this contribution. 


2. Survey of Helios observations in the inner heliosphere 

The results in this section represent a selection out of a more comprehen- 
sive data set obtained by the Helios solar probes during the phase about solar 
activity minimum from late 1974 until the first month of 1976. A complete dis- 
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cussion of the present topics can be found in a series of papers by Rosenbauer 
et al. (1977) and Marsch et al. (1981, 1982a, b, 1983) where a full description 
of the data evaluation and analysis procedure is given as well as an outline of 
the Helios plasma instruments and measurement techniques. The magnetic field 
data used in this paper have been obtained by the TU Braunschweig magnetometer 
experiment (Neubauer et al. 1977). 

The Figure 1 illustrates the variety of shapes and range of functional forms 
of observed proton velocity distributions encountered under various plasma con- 
ditions between 1 AU (top row) and 0.3 AU (bottom row). The three colums pertain 
to low (v < 400 km/8), intermediate speed wind (400 - 600 km/s) and fast solar 
w T ind streSms (v > 600 km/s). The isodensity contour plots represent cuts through 
the distribution in a plane defined by the magnetic field vector (dashed line, 
axis of gyrotropy) and bulk flow direction (VX-axis). Spacing of contours corre- 
sponds to 80, 60, 40, 20 % of the maximum phase space density which defines the 
origin of velocity space. The "core" part of the distributions may be defined by 
this continuous contour line system. Broken contours are spaced logarithmically 
at 0.1, 0.03, 0.01, 0.003, and 0.001 times the maximum. The last contour line 
roughly delineates the instrument one-count level. 

Apparently, solar wind proton distributions come in various shapes and as a 
rule deviate considerably from a Maxwellian. Most of the non-thermal features 
shown here have previously been observed at 1 AU by many spacecraft. Field-aligned 
bulges constituting a total temperature anisotropy T „ > T ± and a heat flux 
appear to be ubiquitous (see also Hundhausen et al. p 1967, p Feldman et al. 1973a) 
with the exception at magnetic sector boundaries (distributions A and G in Figure 
1) where fairly isotropic distributions usually occur. Frequently, the heat flux 
tail attains the shape of a second resolved peak giving rise to a double humped 
ion distribution (see also Feldman et al. 1973b, 1974, Goodrich and Lazarus 1976, 
and Belcher et al. 1981). Double streams are temporarily also observed in the 
alpha particle distributions (Feldman et al. 1973b, Asbridge et al. 1974, and 
Marsch et al. 1982b). A statistical analysis of the speed difference between two 
resolved proton peaks revealed a close correlation of this drift with the Alfven 
speed (Asbridge et al. 1974, 1976 and Marsch et al. 1982a, b) . This observation 
suggests a local regulation of the differential speed by wave-particle inter- 
actions. 

Inspection of the right column shows that in the core part of fast solar 
wind velocity distributions the temperature is larger perpendicular than parallel 
to the magnetic field. This finding may indicate the effects of transverse wave 
heating (Feldman et al. 1973a, and Bame et al. 1975) and of cyclotron resonance 
(Marsch et al. 1982a). Further evidence for local interplanetary ion heating in 
the Helios data is discussed below. In fast solar wind the protons show T A > T ,, 
in the core whereas for alphas one finds T x < T „ in the central part p p 
of the distributions (Marsch et al. 1982a)? This result suggests that the pro- 
cesses shaping ion distributions affect the various ion species in a very differ- 
ent manner. We also recall that in the body of high speed streams alphas move 
faster than protons by about the local Alfven speed (Asbridge et al. 1976, Marsch 
et al. 1982b). In summary, Helios observations yield an entire "zoo" of distri- 
bution functions exhibiting various non— thermal features that become most pro- 
nounced at perihelion (0.3 AU). The association of these diversely shaped distri- 
butions with macroscopic stream structures has been investigated by Feldman et al. 
(1974) and in the paper by Marsch et al. (1982a), that also contains a more 
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Figure 1 Helios 2 proton velocity distributions as measured for various solar 
wind speeds (increasing from the left to the right-hand side) and 
various radial distances (decreases from above to below) . The cuts 
through the three-dimensional distributions are provided in a plane 
defined by the bulk velocity vector (VX-axis) and the magnetic field 
vector (dashed line). Contour lines correspond to fractions 0.8, 0.6, 
0.4, 0.2 of the maximum phase space density (continuous lines) and 
logarithmically spaced to fractions 0.1, 0.032, 0,01, 0.0032, 0.001 
(dashed lines), respectively. The origin of velocity space is defined 
by the velocity of the maximum phase space density and scales are 
given in km/s. 
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detailed discussion of double humped ion distributions and tables of relevant ion 
parameters. Here, we don’t have space enough to discuss all these features. There- 
fore, in the remainder of this section we shall concentrate on the heliocentric 
radial temperature profiles, the so-called adiabatic invariants, and the proton 
heat flux density. 

Before embarking on the discussion of actual data we briefly recapitulate 
the energy equations which are obtained for each ionic species within the frame- 
work of fluid models (Chew et al. 1956) 


si x 

Dt T j" 


3T.„ 


+ 2T.„ b*(V lt Vj) 
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Here the Boltzmann constant is k 
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vective derivative, and V,, = b^\7 with the ^ magnetic field unit vec- 

tor b * B/|b|. The two components of the heat flux tensor are q.n,x defined by 
third order moments. If the random velocity^in the species rest^ frame with velo- 
city V. is denoted by w, one has q.„ = < w„ > n.m^ and q. ± » < w ll w 1 2 /2 > n.m. 

^ normalized dis- J 


where ^ brackets indicate averages J over the 


J J 


tribution function. The field aligned heat flux vector is given by Q, = b (q . «• + 

2q. x )/2. The right hand sides of equation (la, b) comprise the heat ^sources^due 
to^the divergence of the heat flux tensor and also contain heating or cooling 
rates related to Coulomb collisions or wave-particle interactions (index c, w). 
Without these dissipative terms one can combine the two equations with the result 

tjJ 2 2 ..... 

— (T.ii T., / n. ) * 0. By exploiting the frozen-in-field condition and the 

Dt J J 1 - J 

continuity equation we can also write the double-adiabatic equations of state in 
the usual form 


T. l / B = const 
2 

Tj„ (B/n^) = const 


(2a) 

(2b) 


The validity of these equations can be tested by analyzing radial profiles of ion 

temperatures and densities, and the gradient of the independently measured mag- 
netic field strength along individual stream lines. The practical difficulties^ 
encountered in an experimental test of (2a, b) and the assumptions to be made in 
that analysis are extensively addressed in the paper by Marsch et al. (1983). Our 
following discussion is based on heliocentric gradients of ion plasma parameters 
that have been obtained by sorting the data according to the wind speed and by 
averaging over radial distance bins of 0.1 AU width. 


The perpendicular proton temperature T ± is shown versus radial distance 
from the sun in Figure 2 in a double logarithmic plot for various wind velocity 
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classes indicated at the curves, respectively. Within a velocity class, the tem- 
perature increases with decreasing solar distance. At a fixed radial distance, 

T x increases with increasing wind speedy (Burlaga and Ogilvie, 1973). Note the 
P very high temperatures of almost 10 °K at perihelion. Steepest gradignts 
are observed in fast streams. Least squares-fits yield power laws T x ^ r with 
an index “a" ranging between 0.8 and 1.2. These average radial -2 pr °^^ eS are 
flatter than expected for an adiabatic expansion that implies T jyr for simple 
spherical geometry. 


This result becomes more apparent by plotting the proton magnetic moment. 

As can be seen in Figure 3, T ± /B increases with increasing heliocentric distance 
demonstrating that the magnet?c moment is not conserved and that adiabatic in- 
variance is violated. Similar curves are obtained for alphas, however in fast 
streams (v > 600 km/s) the curves are almost horizontal. For protons least- 
squares a fits yield p 'v r a with 0.6 a < 0.9. This trend would be even more 
pronounced if the p magnetic moment were solely based on the “core" distri- 

bution. We interpret Figure 3 as evidence that proton distributions are trans- 
versely heated in the interplanetary medium between 0.3 and 1 AU (see also Bame 
et al. 1975). A more complete discussion of heliocentric temperature profiles and 
adiabatic invariants of solar wind ions can be found in Marsch et al. (1982a, b, 
1983). It should be mentioned that a recent analysis of Voyager plasma data at 
large heliocentric distances of several AU also provided evidence for a ncn-adi- 
abatic evolution of the proton temperature and the occurence of interplanetary 
heating (Gazis and Lazarus, 1982). 

In Figure 4 the second adiabatic invariant of the alphas is shown versus 2 
radial distance in the same format as before. Apparently, on the average T^,,(B/n^) 
decreases during the solar wind radial expansion. Despite large scatter in the 
data caused by large absolute uncertainties in the measured n^ and by spatial in- 
homogeneities and temporal variations the radial trend in Figure 4 is believed to 
be statistically significant. This finding may be interpreted as evidence for a 
non-adiabatic evolution of T ft , the alpha-particle parallel temperature. T^n 
seems to decline more rapidly than expected for adiabatic cooling. This result 
could well be related to the radial evolution of the ion differential velocity^ 
which is observed to trace the Alfven velocity (that roughly scales like ^ r ) 
and which thus decreases during the radial expansion. 

In equation (la, b) the terms giving rise to a non— adiabatic behaviour of ion 
temperatures are the divergence of the heat flux tensor and Coulomb or anomalous 
(wave-particle interactions) collisional transfer rates for T.„ , A . For Coulomb 
collisions analytic expressions for U temperature exchange” between the ions 
can be found in the classical article by Braginskii (1965) on collision dominated 
transport. Observations in slow, cold solar wind indicate the possibility of an 
effective coupling between the various kinetic degrees of freedom for protons and 
alphas by Coulomb collisions (Feldman et al. 1974b, Neugebauer 1976, Neugebauer 
and Feldman 1979, Griinwaldt and Rosenbauer 1978, Neugebauer 1981, Marsch et al. 
1981, 1982b). High speed wind data on the other hand indicate, though, that Cou- 
lomb collisions are much less important in hot, fast solar wind streams. Without 
invoking wave-particle interactions under these conditions there only remain the 
heat flux terms as a possible source of ion thermal energy. 

The Figure 5 displays the proton heat flux density Q = (q „ +2q ,)/2 
versus radial solar distance in a semi- logarithmic plot for various solar 
wind velocities. Note that the variations in extend over two orders of mag~ 


359 









nitude. Between 0.3 and t AU Q ranges from several 10~ 4 to about 10~ 2 ergs/cm 2 s 
Least squares fit yield Q ~ p r with 3.8 < a < 4.7 and Q ^3-4 10 -4 ergs/cm 2 s 
at 1 AU for intermediate p and high speed solar wind. Alpha p % particle heat 
flux densities turn out to be almost an order of magnitude smaller (Marsch et al. 
1983) which means that q^n, x can be entirely ignored in equation (la, b). For 
protons one finds that q . is much smaller than q „ (implying Q ^ q „) . Thus in 
equation (tb) the terms with q . can also be ^ neglected. In^View^of these 
results, one needs wave-particle p interactions in order to account for the radial 
course of the proton magnetic moment in Figure 3. This conclusion seems to be 
compelling if one considers the temperature anisotropy in the high speed proton 
distributions of Figure 1. An anisotropic local heat source 3T. x /3t | in terms 
of wave heating should be invoked to explain this detail in J proton w distributions 
and to further explain the differences in the shapes of the alpha and proton 
distribution functions (Marsch et al. 1982a,b). 

Concerning the parallel temperature T.„ one can roughly estimate the strength 
of the interplanetary heat source represented by — V*(q.,,b). By integrating this 
quantity over the radial heat flux profiles given by the ^ least squares fits in 
Figure 5, Marsch et al. (1983) found that the temperature increase resulting from 
a degradation of the heat flux cannot account for the non-adiabatic radial evo- 
lution of the proton temperatures. For the alphas the result of Figure 4 even 
suggests the need for an interplanetary heat "sink" or effective redistribution 
of thermal energy in order to explain the stronger than adiabatic cooling of T „ 
between 0.3 and 1 AU. a 

We conclude with a presentation of Helios 2 data that pertain to the phase 
of the solar cycle inclining to solar maximum. In contrast to the typical solar 
minimum plasma conditions, during a few hours on day 117 in May 1978 in the peri- 
helion at 0.3 slow solar wind was observed associated with pronounced ion differ- 
ential speeds (AV ^ 100 km/s) and temperature anisotropies (T x > T „) and with 
unusually intense F Alfvenic type wave activity (Marsch et al?1981)? Figure 6 
shows from top to bottom proton speed, ion differential speed (Alfven speed is 
indicated by points) , the dot product of AV and B , and the proton and alpha 
temperature anisotropies T.„/T. x . The core -<Xp temperature ratio is marked by 
points. Note the low wind ^ ^velocity. It stays low' at about 400 km/s despite 
large Alfvenic-type velocity fluctuations. Simultaneously, the alphas move faster 
than the protons by more than 100 km/s. They move at about the local Alfven speed 
along the interplanetary magnetic field (panel three). In the core part of the 
proton distribution one finds T x > T „ whereas for the alphas the reversed situ- 
ation occurs. Despite large P P scatter in the data due to worse counting 
statistics, on the average T qA < T a „ is observed. These contrary features have 
earlier been found to be characteristic of distributions in fast recurrent streams 
as discussed in Figure 1 . In our opinion, this result is a key observation towards 
understanding the mechanism of solar wind heating and acceleration. Obviously, 
these ion distributions exhibited properties that have been considered typical of 
high speed wind originating from coronal holes during solar activity minimum. How- 
ever, the unusual plasma investigated here was embedded in "normal", low speed 
wind. The ambient plasma showed the typical signatures of slow wind as discussed 
in the beginning of this section. 


361 


800 


Figure 6 




RCA.U. )- O.290 0.200 0.200 0.200 0.200 0.200 0.200 

LAT °- 1.02 1,00 1.00 1.01 1.04 1.00 1.00 

LONG 0 - 103.1 102.0 100.7 100.8 100.2 100.0 101.0 

DATE-27.4.78, DAY 117 

HELIOS 2 



Six hours of Helios 2 ion data 
from day 117 in 1978, measured 
at 0.29 AU. Reading from above, 
the proton bulk speed, the ion 
differential speed (and indi- 
cated by points in the same 
panel, the local Alfven speed), 
the cosine of the spatial angle 
between Av and the magnetic 
field B^, ~ and finally, in 
the lowest panels, the proton 
and helium ion total tempera- 
ture anisotropies are shown 
(the core anisotropy is marked 
by points) . 


3. Theoretical models and discussion 

Certainly, the. radial evolution of the detailed solar wind ion distributions 
cannot be explained within the framework of magnetohydrodynamic fluid theories 
but rather requires a kinetic treatment. Such an approach to the problem has yet 
not been tried. However, a vast amount of literature exists concerned with ana- 
lysing the stability of distributions observed in the interplanetary medium. A 
comprehensive review on microinstabilities is given by Schwartz (1980) and on 
various kinetic aspects by Feldman (1979 a, b) . A theoretical discussion of wave- 
particle interactions with emphasis on MHD-waves is also included in Barnes' 
(1979) review on wave turbulence in the solar wind. 


362 



„ In instability calculations particle distributions are shown to relax 
locally to a state that corresponds to less free energy than the original state. 
Reordering of internal ion energy is accompanied by excitation of various wave 
modes where growth time and dispersion characteristics sensitively depend on the 
shape of the actual distribution. As demonstrated in Figure T, the main free 
energy sources for local wave excitation are represented by temperature aniso- 
tropies, heat fluxes and double streams, and also by the differential proton 
alpha-particle motion that was discussed in some detail before. Instability cal- 
culations have been performed for given this free energy in non-thermal ion 
distributions. However, no attempt has yet been made to explain, e.g. double ion 
streams in the first place by dynamic processes in the corona or during the 
initial phase of solar wind expansion. Feldman et al. (t974) proposed conceptual 
models in terms of double ion streams originating from time variations in the 
corona or from the interpenetration of fast into slow ambient plasma. However, 
these ideas are highly speculative and at present very little is understood con- 
cerning the origin of non— thermal features of ion distributions. 

Goodrich (1981) investigated the kinetic effects of Alfven wave pressure on 
the ion distributions. By solving the. relevant diffusion equation he could show 
that the kinetic nature of Alfven wave— particle interaction is reflected in 
distortions of the distribution function that give rise to well known fluid 
results like bulk acceleration. However, these model calculations are far from 
accounting for the detailed shape of the observed distributions. Modeling the 
radial evolution of ion distributions in terms of convected Bi-Maxwell ians re- 
presents an intermediate step between a fluid and fully kinetic treatment. 

Namely, by calculating the temperatures and speeds relative to the solar wind 
frame one readily knows the structure of the distribution functions that are 
entirely specified in terms of their moments. Such model assumptions have almost 
universally been adopted for plasma instability calculations. 

Along these lines Marsch et al. (1982c) have developed a self consistent 
fluid-type model in order to explain the ion temperature anisotropies and differ- 
ential speeds by means of wave heating and acceleration of solar wind ions by 
cyclotron resonance with ion-cyclotron and magnetosonic waves propagating away 
from the sun along the interplanetary magnetic field. Local temperature and 
momentum transfer terms based on the quasilinear expressions for wave-particle 
interactions in weak electromagnetic turbulence (see also Dusenbery and Hollweg, 
1981, and McKenzie and Marsch, 1982) have been employed in the fluid equations. It 
was shown by Marsch et al. (1982c) that the magnetosonic waves are capable of 
accelerating the minor ions to differential speeds of the order of the local 
Alfven speed. Furthermore, it was demonstrated that interaction with ion-cyclo- 
tron waves can lead to an equalization of ion thermal velocities and the observed 
anisotropy of proton temperatures in close qualitative accord with the obser- 
vations. 

However, the required wave intensity exceeds the values observed (at 0.3 AU) 
by at least one order of magnitude, which is a serious problem all theories are 
faced with that rely on resonant wave-particle interactions at about the ion gyro- 
frequencies (Schwartz et al., 1981). This problem might be circumvented by invoking 
a cascade process from the low frequency Alfven waves to the resonant frequency 
regime (Isenberg and Hollweg, 1982) in order to replenish the intensity at higher 
frequencies. Still, the physics of tapping the energy reservoir of Alfven waves 
by a cascade remains rather obscured. In addition, the assumption of rigid model 
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distributions in the quasilinear expressions represents another problem since the 
plasma may in reality not conform to these model constraints simply by changing 
the shape of the distribution function. This effect can be crucial in case the . 
interactions take place in the tails of the distribution. Using real distribution 
functions measured on Helios, Dum et al. (1980) showed that for the ion acoustic 
mode and whistlers the growth and damping characteristics depend rather meticu- 
lously on the detailed shape of the particle distributions. At present a fully 
kinetic treatment of the problem is still lacking, although it. appears to be 
needed for an explanation of the radial evolution of ion distribution functions. 


Finally, we briefly mention an alternative explanation of the large perpen- 
dicular ion temperatures in terms of Landau damping of lower hybrid waves. Marsch 
and Chang ( 1983) have demonstrated that the broad band low frequency electro- 
static noise frequently observed in the disturbed solar wind at interplanetary 
shocks (Coroniti et al., 1982; Kennel et al., 1982) and in high speed sola * 
streams (Beinroth and Neubauer, 1981, and Gurnett et al., 1979) may generaliy hav 
a lower hybrid component. These modes accompanied by hybrid-like whistler 
can be excited by resonant halo electrons in the heat flux tail of the e ectron 
distribution (Rosenbauer et al., 1977; Feldman et al., 1981). Model calculations 
based on convected Bi-Maxwellian electron distributions (Marsch and Chang, 1^-3) 
have shown that electromagnetic lower hybrid waves at several - v ' e i' , 

can energize solar wind ions transverse to the magnetic field, . since . 

waves propagate almost perpendicularly. By means of Landau damping of lower hybrid 
waves the proton distributions could thus attain large perpendicular temperatures 
in the core. Furthermore, the effective coupling between protons and supratherma 
electrons, being simultaneously in Landau resonance with the waves, should be of 
major importance for understanding solar wind transport and the regulation ot tne 
electron heat flux. 


A. Conclusions 

In this paper we have discussed some selected aspects of the. radial evolution 
of ion velocity distributions that exhibited collisional and collision- ess e 
haviour. In slow, cold, and dense solar wind there is evidence for a collisional 
redistribution of internal ion energy (see also review by Neugebauer, 1981). n 
contrast, high speed ion distribution indicate significant influence of waves . 
that render the solar wind expansion non-adiabatic. Various processes shaping. ion 
distributions by means of microinstabilities and models that attempt to describe 
the radial evolution of internal ion energy have been proposed in the literature. 
Two basic lines, that a theoretical description of the solar wind expansion 
principle could follow, are sketched in Figure 7. The standard fluid approach has 
widely been used in many theoretical papers at this conference (see. also the SWA 
review by Cuperman, 1981). However, very little work has been done in terms of 
kinetic equations (Griffel and Davis, 1969, and Eviatar and Schulz, 1970) that^are 
solely appropiate to describe the radial evolution of. detailed ion 
and internal energy and to take dissipation properly into account. At presen , 
observational knowledge is far ahead of a coherent theoretical understanding of 
solar wind microprocesses and transport. Therefore, a kinetic descrip 1 ‘ , 

wind ions that incorporates collisions and wave-particle interactions is urgently 

required . 
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Figure 7 Schematic sketch of two possible basic approaches to a theoretical 
description of the solar wind expansion. 
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